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Abstract: We present a new model of BEC dynamics based 
on strong confinement near the ground state. The model pre- 
dicts oscillations in a two-component condensate, based on 
interference of non-spreading wave packets moving within a 
pair of tilted nearly square potentials. The oscillations are 
similar to those recently reported for a magnetically trapped 
*^Rb condensate, and the model's predictions give good quan- 
titative agreement with the experiments. 
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Since its first application to quantum Rabi oscillations 
in the two- level cavity QED context, collapse and re- 
vival analysis |Q] has been used for the understanding of 
a widening class of recurrent physical phenomena. The 
timing and nature of revivals have been employed in in- 
terpreting the dynamics of a variety of quantum systems 
1^,^ , but a number of cases provide puzzling features that 
are considered open questions 

Recent elegant experimental observations of Rabi-type 
oscillations in a two-component Bose-Einstein conden- 
sate of ^^Rb 1^ offer a new example that has resisted 
analysis. Multi-component EEC's |^,^ in the presence 
of external driving fields have provided a rich spectrum 
of physical responses to challenge current theoretical un- 
derstanding. In the two-component case of interest here 
the theoretical situation has been described by var- 
ious groups P,p^. In particular, differential shifting of 
the trap centers for two species has led to interesting the- 
oretical predictions such as a "condensate dressed state" 
[ pl| , Josephson-type oscillations and periodic modulation 
of Rabi oscillations , all of which originate in analogies 
with two-level quantum optical systems. Despite heuris- 
tic arguments using "twists" of an SU (2) order parameter 

, it is fair to say that the fundamental dynamics of the 
observed oscillations are not understood. 

In this Letter, we go beyond revival analysis to present 
a new quantum mechanical picture of the physics based 
on a strong confinement model, which explains ex- 
perimentally observed features quite well. We show 
that the coupled Gross-Pitaevski (GP) equations for 
the two-component Bose-Einstein condensate can be de- 
coupled by a time-dependent dressing transformation. 
The transformed Schrodinger equations describe local- 



ized wavepackets that accelerate in opposite directions 
and bounce back and interfere after hitting the walls of 
the confinement. We find that the spatial and tempo- 
ral behavior of the two BEC's are determined by these 
wavepackets and the experimental observations that have 
been reported can be explained in detail in terms of their 
interference. Our model gives quantitative predictions 
that also agree well with numerical simulation of the cou- 
pled GP equations. 

The two components of the condensate, "01 ^-nd "02, 
obey GP equations that are coupled to each other via 
the "drive" parameter f2, which plays a role similar to 
the Rabi frequency of a driven two-level system: 

= Hoi'i + t^{ri\iPi\^ + \lP2\^)lPl + ^Vi + (1) 

«^ = ^^0-02 + fJ-il-tpll"^ + r2\lp2\'^)'>Jj2 - ^V'2 + ^-01, 

(2) 

where Hq = —\-§^ + is the bare Hamiltonian for the 
harmonic trap. For simplicity we present here only a one 
dimensional model. We have used the axial oscillation 
frequency Uz to express units: Tiujz is our unit of energy 
and l/uJz and dz = ^/'hj mtOz are the units of time and 
distance. The coefhcient ^ = 4:'Kai2N/dz measures the 
strength of interaction between two species of BEC where 
N is the total particle number and ai2 is the interspecies 
scattering length. The small difference reported |^ for 
the ratios in ri — ai/ai2 ~ 1.03 and r2 = a2/ai2 ~ 

0.97, among scattering lengths has little effect on the 
subsequent discussion so we assume that ri = r2 = 1 in 
the following. 

The detuning S and the interspecies coupling Q, both 
measured in the unit lUz, can depend on the axial posi- 
tion z. We will make assumptions of no detuning {S = 0) 
and a linear z-dependence of coupling: O = fio + 2/9z, 
consistent with experiment |^] . The nonvanishing detun- 
ing case will be discussed later. Now we introduce the 
time-dependent dressing transformation 

0± ^ (V-i ± 02)e±^""*/^ (3) 

which gives two decoupled GP equations in a slow time 
frame: 
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where the potentials V± are given by 

F±(z) = lz2^^(|0+P + |0_|2)±/?z. (5) 

The first crucial observation is that although these po- 
tentials are spatially very nonlinear, they are nearly in- 
dependent of time if the system is started near its ground 
state. This is due to the strong and relatively tight 
nonharmonic confinement provided in the present case 
by the GP interaction. In principle, the potentials de- 
pend on time through the total density of the conden- 
sate, ptot = IV'iP + IV'2p = ^(I'/'+P + I'/'-P), but ptot 
changes very little during hundreds of Rabi cycles that 
are well-resolved in the experiments. This is evident in 
the relatively very rigid spatial shape of the total conden- 
sate density seen in the series of snapshots in Fig. 1. To 
obtain the snapshots we solved the coupled GP equations 
using the combined split-step Fourier technique familiar 
from fiber optics theory (where the GP equation is called 
the non-linear Schrodinger equation), modified to include 
two components |l3| 

The second crucial observation, also evident in Fig. 1 
and noted previously |^ , is that ptot maintains a Thomas- 
Fermi-like spatial distribution. We have found that this 
allows the confining forces on the packets to be converted 
into boundary conditions, leaving relatively free (f>± pack- 
ets to accelerate slowly in opposite directions at the rate 
ztf3 along the bottoms of two potential wells that are 
nearly square but tipped in opposite directions, as shown 
in Fig. 2. We will now demonstrate that this simpli- 
fied picture is powerfully effective in unraveling relatively 
complicated condensate evolution. 

If, as in the JILA experiments, a condensate at nearly 
zero temperature is prepared initially in the pure 1 com- 
ponent, i.e. -i/jilzjO) = 'ipgiz), ip2{z,0) = 0, we have 
equally distributed ip's : 0+(z,O) ~ 0_(z,O) = ipg{z). 
As these ± probability packets fall along the bottom of 
the well toward opposite walls they must interfere since 
(/)+ and (j)- are coherent parts of the same two-component 
condensate. An effective interaction between the decou- 
pled ± packets comes from their interference to produce 
the observed density It/^iI^: 

«^IV'sP + ^5R(e»"oV;0-), (6) 

where 5ft denotes the real part. 

First we show additional results of direct numerical 
integration of the original two nonlinear GP equations 
(^. The calculated temporal evolution of the space- 
integrated population of component 1 is shown in Fig. 
3. The corresponding experimental data is given in Fig. 



4-a. Rapid oscillations at the drive frequency Q.q are ev- 
ident, as the population cycles between components 1 
and 2. An apparent collapse in region A, and a revival 
or perhaps a pair of them, is seen in regions D and E. 
Small features at B and C are also noted. Now we will 
use the strong confinement model to obtain a detailed 
understanding of the curve. As we will see, revivals in 
the conventional sense play no role. The correct expla- 
nation is in some respects simpler, and capable of wider 
applicability as we mention at the end. 

Now we return to our simplified picture, according to 
which the main slow motion is acceleration with each 
packet moving toward an opposite wall. This main slow 
motion can also be factored out of the quantum ampli- 
tudes by defining 

X±{z, t) = (t)±(z, t) exp{±i(3zt + i(3^ty6), (7) 

and changing the coordinates into those of the two 
accelerating frames: T — t, Z± = z ± pt^ /2. 
The packets x± satisfy the free Schrodinger equation, 
{id/dT+\d'^/dZ'^)x± = 0, with initial conditions 
X±{T — Q) = ipg, and so can be integrated directly to 
yield 

Since the initial tpg is very smooth, this shows that the 
two x± move apart with a combined quantum mechanical 
dispersion which is relatively small in our Thomas-Fermi 
like configuration. The interference term in Eq. (^) is now 
given by ^Sft (exp[zrioi + 2i(3zt\x%X~') with x± as given 
above. 

Thus we predict spatial oscillations with wavelength 
Lq = TT / (3t to develop in the population of condensate 1, 
dramatically evident in Fig. 1-B. The interference occurs 
in the region where the ± packets continue to overlap, a 
shrinking region of length L ^2R- I3t^, where the fuU 
axial length is 2R. These spatial interference oscillations 
have nodes that can be monitored experimentally [see 
Fig. 4-b] and the node formation has been called "twist- 
ing." Various effects are predicted and labeled in Fig. 3 
as follows. 

A. Collapse. The number of nodes in existence at time 
t is fixed at n = L/Lq = /3t(2i? - /3i^)/7r. By about the 
time of the second twist there is a good balance of the 
two components in the condensate (e.g., ~ l'02p ~ 
0.5), which shows up in the data as the completion of 
the "collapse" at the end of region A in Fig. 3. From 
n = 2, with the parameters used in the calculation [R = 
8.53, j3 — 0.25), one easily finds for the collapse time 
tc ~ 1.5, in excellent agreement with Fig. 3. 

B. Twist resolvability. Our formula also predicts a 
maximum twist number hm, and says that it occurs at 
time tM = •\/2i?/3/3, after which the ± packets are falling 
too fast to have resolvable spatial interference. This gives 



2 



Um = ^/iWWp^^. or tM = 4.77 and um = 4.3, also 
in good agreement with the observed behavior in region 
B. 

C. Bounce time, t ^ ts- When the centers of (j)± reach 
the walls at time tB — y/2R/f3 ( = 8.26 in our case), 
the (f)± become maximally separated and have a finite 
but minimally overlapping interference region. Eq. (j^) 
shows that in the near- wall regions, the densities jV'iP 
and \tp2\^ reduce to about half of the initial ground state 
density. This cannot be seen in Fig. 3, but is confirmed 
by examining snapshots of the spatial behavior of the 
two condensate components. Five such snapshots were 
already given in Fig. 1, and Fig. 1-C confirms that 
at the time of the wall bounce, in the regions near the 
walls, the two components are not interfering but each 
contain about half of the total density. Of course, due 
to reflection at the wall and subsequent self-interference, 
there arise complicated motions for each A more 
complete analytic description of the interference between 
(f>± can be made by solving the Schrodinger equation (Q) 
directly in terms of Airy functions. 

D. Packet recovery, t = tR. After bouncing from 
the strong confinement, the packets 0± return approx- 
imately to their initial positions. As they return, os- 
cillations in the overlapping region spread out over the 
half-ground-state density region with increasing ampli- 
tude and wavelength. When the number of nodes reaches 
three, with one central peak and two side peaks, oscil- 
lations in the half-ground-state region interfere most de- 
structively, thereby minimizing the population of compo- 
nent 1 (Fig. 1-D). This happens before completing one 
round trip by the amount of time required to develop 
three nodes in the initial stage. However, the external 
coupling drive quickly switches back from the minimum 
population to the maximum population for condensate 
1. Thus, the maximal recovery occurs at = 2<b — At 
where Ai is obtained from (3At{2R - (3{Mf)/iT = 3. 

In our case. At — 2.4 and t^ — 14.1. This shows 
that the maximum revival of population docs not come 
with the original shape. In fact. Fig. 1-D shows that it 
rather has the shape of a dark soliton. This is also in 
agreement with the reduced contrast of about 60 percent 
of the original. 

E. Completing one cycle, t — tp. At the classical re- 
turn time, tp = 2tB^ the recovery of population of con- 
densate 1 happens but in a smaller scale than in case D, 
as shown in Fig. 1-E. 

In closing, let us return to an earlier point. In the 
spin-boson language of revival analysis, a vanishing de- 
tuning [5 ~ Q in the evolution equation (^) , for example) 
corresponds to exactly zero Zeeman splitting of the two 
spin states, removing any possibility of resonant coinci- 
dence with the boson (oscillator) transition frequencies. 
This is why revival anaysis is misapplied to the existing 
data. However, although the BEG oscillator is very an- 
harmonic under its strong confinement, one suspects that 



it could support a near-resonance with a non-zero 5, if i5 
were large enough. In such a case one would have an in- 
triguing combination of bouncing packets and quantum 
revivals (which themselves are of course known to have 
packet interpretations). Since the energy spacings in the 
confinement are small, 5 would not have to be very large. 

In this Letter, we have introduced a quantum mechan- 
ical packet model which provides an excellent qualitative 
description of the previously unexplained temporal be- 
havior of oscillations in two-component driven-BEC ex- 
periments. One clearly sees that the strong-confinement 
bouncing packet model permits easy analysis that pre- 
dicts all five distinct events in the temporal record as 
well as predicting their appropriate time scales. Figs. 1- 
3 compare these events with numerical simulations and 
with the experimental record in Fig. 4. 

It is clear that this model is not restricted to the two- 
component case, and it will be interesting to extend our 
work to the three-component BEG recently reported 
and partially analysed |p^ . 
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FIG. 2. Potentials V±. Since the total density maintains 
the Thomas-Fermi shape, the potentials have a flat bottom 
while the wall is the upper portion of a harmonic potential. 




FIG. 3. Plot of the population of condensate 1. The 
thick line represents the case with no constant coupling drive 
{Qo ~ 0) while the fast oscillating thin line represents the 
case with flo = 28.85. Five distinct events are specified at 
the calculated times, which are described separately in the 
text. 
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FIG. 1. Calculated density profiles of condensates at five 
distinct times that are the keynote instants during evolution 
noted in Fig. 3. Solid and dotted lines represent condensates 
1 and 2, and the dashed line represents the total density ptot- 
Note the strong time-independence of ptot- 
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FIG. 4. Experimental data from ref. [5]. The modest dis- 
crepancy in time registry copared to Fig. 3 is due to our 
simplified ID modeling and also to an experimental depar- 
ture from 5 = 0, evident because the mid-line in part 4(a) 
shows a component 1 fraction lower than 50%. 
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